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Abstract 



We calculate the baryon asymmetry generated at the electroweak phase transition in the 
minimal supersymmetric standard model, treating the particles in a WKB approximation 
in the bubble wall background. A set of diffusion equations for the particle species relevant 
to baryon generation, including source terms arising from the CP violation associated with 
the complex phase S of the /i parameter, are derived from Boltzmann equations, and solved. 
The conclusion is that 5 must be ^,0.1 to generate a baryon asymmetry consistent with 
nucleosynthesis. We compare our results to several other recent computations of the effect, 
arguing that some are overestimates. 



Although theories to explain the baryon asymmetry of the Universe (BAU) abound, only 
electroweak baryogenesis has the hope of short-term testability. To calculate the asymmetry 
in detail one needs a specific model of electroweak physics beyond the standard model, and 
supersymmetry is the most popular idea for such new physics. Several groups have recently 
estimated the BAU coming from the minimal supersymmetric standard model (MSSM) at 
the electroweak phase transition (EWPT), assuming it is first order and so proceeds by 
bubble nucleation |y-||. Despite the fact that the physics of the mechanism is essentially 
agreed upon — particles in the plasma interact with the bubble wall in a CP- violating manner, 
leading to a chiral asymmetry of quarks in front of the wall, which in turn biases sphalerons 
to produce the BAU — there is considerable disagreement on both methods and results. 

While it is straightforward to formulate diffusion equations for particle transport away 
from the wall, it is less obvious how to derive the CP-violating source terms (or equivalently 
the boundary conditions at the wall) for these equations, and it is here that most of the 
difficulty lies. Several methods have been proposed in which the source terms are first 
derived and then inserted into a set of diffusion equations, but without any first principles 
justification for the prescription used. In the present work we generalize a method which was 
introduced for the treatment of the non-supersymmetric two doublet model in 0], in which 
the diffusion equations and the source terms are derived together starting from a general set 
of classical Boltzmann equations. This provides a systematic and controlled approximation 
valid for bubble walls significantly thicker than the inverse temperature, which is the case 
in the MSSM. Here we will describe the results as concisely as possible; finer details will be 
presented in a longer publication ||. 

The method is based on a WKB quasi-particle approximation to the particle dynamics 
and interactions in the background of the bubble wall, which is an expansion in derivatives of 
the Higgs field VEV's that vary continuously in the wall. Such an expansion is appropriate 
because it well-describes the majority of particles in the plasma, whose wavelength (~ 1/T) 
is much shorter than the thickness of the wall (~ 20/T). When there is CP violation 
in the theory, such as is provided by the complex [i and A t parameters of the MSSM, 
this approximation shows that particles and anti-particles experience a different force when 
interacting with the bubble wall, both dynamically as free WKB particles and in their 
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interactions with other particles in the plasma. These two effects — the classical dynamical 
force on the particles, and the CP biasing of their collisions — appear in the Boltzmann 
equations for the WKB particles in a quite straightforward way, through the force term 
and collision term, respectively. These equations determine the distribution of left-handed 
quark number (q^) in front of the bubble wall, which is the source term for baryon number 
generation. An appropriate fluid-type truncation leads to a set of diffusion equations of a 
familiar form, together with the source terms, which predict the qi and hence the baryon 
asymmetry. 

In what follows we will concentrate on charginos as the source of the chiral quark asym- 
metry, because we find they make the dominant contribution. The corresponding analysis 
for the subdominant squark and quark contributions will be presented in ref. ||. To find 
the effect of CP violation on charginos at the bubble wall, we first solve the Dirac equation 
in the WKB approximation. 

WKB approximation for charginos. The mass term for charginos can be written 
as ^/rM x ^/ l + h.c, where in the basis of Winos and Higgsinos, ^/r = (W^,hf R ) T , = 
(W^! ^A,l) T the mass matrix is 

with the spatially varying Higgs field VEV's t>j. The corresponding Dirac equation is 
{ilzlodt + id z — 7 3 (M x Pk + M^Pr))^! = 0. In the WKB approximation the local eigen- 
states of mass, energy and helicity are related to the flavor states through a transformation 

* L \ = ( V °)( cosh W e l Mnh(X) \ / ^ \ +i rp(z)dz 
^ R ) \ U J \ e~ ie sinh(X) cosh(X) J \ # T ) ' 1 ' 

Here U and V are the unitary transformations that locally diagonalize the mass matrices, 
WM X V = M x and V^M X U = M* x , and X is defined by tanh(2X) = \M X \/E, which 
determines the mixing between the chirality and the helicity states. 8 is a diagonal matrix, 
consisting of the phases of the eigenvalues Ai x , whose spatial variation is responsible for all 
the CP violating effects we will be concerned with. Finally, p(z) is also a diagonal matrix 
that gives the local momentum of the two mass eigenstates. The explicit expressions for 
the mass eigenvalues, unitary matrices and CP- violating phases can be given in terms of the 



2 



quantities 



m 2 = (ml + |/i| 2 + u\ + ul)/2; Ui=gVi/\^2; 
A = (m 2 , — \fi\ 2 — u 2 + ul)/2; a = m 2 Ui + n*u 2 
A = (ml — |/i| 2 + u\ — ul)/2; a = m 2 u 2 + 

A = (A 2 + |a| 2 ) 1 / 2 = (A 2 + |a| 2 ) 1 / 2 ; (3) 

then |.M Xi ±| 2 = m 2 ± A, and 

U = I ( A + A ~ a ) V = I ( A + A ~ a ) (4) 

V /2A(A + A) I A + A J' ^2A(A + A) I A + A y ' U 

and 

|-^x,±| 2 ^± = ip^clm (/i)m2 cos/3sin/3MC? 2 M/A, (5) 



where 9 2 denotes the spatial derivative normal to the bubble wall, t> c = v^ 2 + f 2 is the 



magnitude of the VEV at the critical temperature, u = yuf + u^ and tan/5 = v 2 /vi. In 
what follows the phase 9 will never be needed except in the combination shown above. 
Finally, to lowest nontrivial order in derivatives and CP-violation, the local momenta are 



p(z) = sign(j> 2 ) {\J E 2 — |-M X | 2 =F sinh 2 (X)<9 2 #) , for particles with helicity ±, and in the 
Lorentz frame where the particle has no momentum parallel to the bubble wall. For the 
antiparticles, the same equation holds but with an overall change of sign for the CP- violating 
part. 

The important thing to notice in these expressions is that the CP-violating angle 9 
appears in two places: in the particle momenta (as a derivative) and in the transformation 
matrices to the WKB basis. The former leads to the dynamical force effect and the latter 
to source terms of the "spontaneous" type, which have previously been discussed in the 
literature in various contexts |5) . In what follows it will be convenient to have the dispersion 
relation for energy in terms of momentum, for the two mass eigenstates v. 

d z 9i 



E± = \jp 2 + \mi\ 2 ±sign(p 2 ) — ( ^p 2 z + \m { \ 2 - \p z \ ) = E ± AE; 

2\/.P 2 + \rrii\ 2 



± = + for R, L, - for L, R. (6) 
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To derive this we used the identity sinh 2 (X) = [\Jp 2 z + \m\ 2 — \p z \j /{2\p z \) and transformed 
to a general Lorentz frame with nonzero momentum parallel to the wall. 

Diffusion equations with sources. Our starting point is the Boltzmann equation 
(d t + x ■ dg + p ■ dp)fi = C[fi] for the local distributions fi(p, x) of the WKB quasi-particles 
with the on-shell dispersion relation just derived. The on-shell approximation is valid since 
the particle widths V are small compared to their energies. To truncate these equations we 
use the ansatz 

for the distribution functions in the rest frame of the plasma, where £7j contains the per- 
turbation AEi ~ 89, i.e., we allow a chemical potential and velocity perturbation in each 
fluid. The validity of the ansatz requires that v w < L/3D @, in order that whatever 
other perturbations are produced are indeed small compared to those of the chemical po- 
tential, the latter being damped only by the slower inelastic processes. The velocity per- 
turbation is required to model the anisotropic response to the force term. Expanding 
to linear order in the perturbation AEi and taking /ij and to be of this order gives 
f±d z ((v w - p z /E )(p z Vi + fii) - v w AE^j = -C [fi], where f± = df±/dE denotes the deriva- 
tive of the unperturbed Fermi-Dirac or Bose-Einstein distribution, f'± = (e^ ± Then 
averaging over momentum, weighting the equation respectively by 1 and p z , we obtain the 
two coupled equations 

-tf/EoM + Vvd = (Ci>; (8) 
v w (p 2 z H-(p 2 z /E ) f i' l -v w ( Pz AE:) = (C tPz ) (9) 

where (•) = / d 3 pf±(-)/ J d 3 pf' + = Ki J d 3 pf±(-)/ J d?pf± and the primes in eqs. denote 
d z . Here we have divided by the distribution function for a massless fermion even if the 
numerator is for a boson, because this ensures that the collision integral on the r.h.s. of 
the Boltzmann equation gives the same value for the rate of a given process, regardless of 
whether one is writing the transport equation for a boson or for a fermion involved in that 
process. The result of this choice is to introduce the factor of Ki on the l.h.s. of the Boltzmann 
equation, which is 1 for fermions and 2 for bosons, in the massless approximation. 
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To evaluate the collision terms in eq. one must substitute the distributions (0) into 

Sdn j 8^\ Vi + V2 ...- Vn „ l - Vn )\M\\f i f 2 ---{i±f n - 1 )^ 

which vanishes for the equilibrium distibutions. Thus the collision terms are linear in the 
perturbations: (Cj) = J2j ^f^j^j and (p z Ci) = (p 2 z ) J2j ^ij^ij^i, where T d ' e are respectively 
the inelastic (decays) and elastic interaction rates for processes coupling states i and j. 
There is further another important contribution to (d) coming from the fact that the 
interactions of the flavor eigenstates, when reexpressed in terms of the local mass eigen- 
states, have factors of e %e ^ coming from eq. (Q). In fact that equation is derived in the 
rest frame of the bubble, so that in the wall frame e t6 ^ — *> e l9 ( z - v ^ t ) ( which to leading 
order in derivatives can be written as e l ( 9 ( z )- v ^ dd ) . The effect is to spoil energy conser- 
vation of the interaction so that 5(Ei + E% . . . — E n -\ — E n ) — > S(Ei — v w d9 + E 2 . . . — 
E n _i — E n ). This gives an additional term in (Cj) ~ v w d9, which is of the type referred 
to as a spontaneous baryogenesis source term. In the case of Higgsinos, the interaction 
yqi3Rh2 can be rewritten as 2/<7l£r(V2+ sinh(X + )e j6,+ \l/| + + V2- sinh(X_)e ie ~\I / |_) leading to 
(Cr ) = v w T y \V2±\ 2 (smh 2 (X±))d9±, where ± is whichever sign is needed to make the as- 
sociated mass eigenstate go smoothly into the pure Higgsino state in front of the wall (see 
below) . 

Finally, the two coupled equations can be reduced to a single one by differentiating (Q) 
and eliminating Vi in favor of //j. Defining the diffusion constant Di = (p 2 z / E) 2 /((pf)Tf) and 
ignoring the ratio of inelastic to elastic scatterings, we obtain the diffusion equation 

-Ki(Din'l + v w fii) + rf^/jij = Sf, 

to leading order in the wall velocity. 

Solution of diffusion equations. From this complicated network of equations coupling 
all species in the plasma we need to determine the distribution of left-handed fermions, 
which bias the unsuppressed baryon-number-violating processes in front of the wall. It can 
be simplified by using conservation laws and neglecting processes which are too slow to play 
any role on the relevant timescales. Parametrically this means comparing the time spent 
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by a particle diffusing in front of the wall ~ D/v^ with the decay time The processes 
we do take account of are then (i) the supergauge interactions, (ii) those described by the 
following terms in the interaction Lagrangian 

V y = yh 2 u R q L + yu R h 2L q L + yu* R h 2L q L - y/j,hiq* L u R + yA t q L h 2 u* R + h.c. (11) 

and (iii) strong sphaleron interactions. Taking the supergauge interactions to be in equi- 
librium, the chemical potentials of all gauginos are zero and the chemical potential of any 
particle is equal to that of its superpartners. Defining a reduced set of chemical potentials for 
each generation and chirality of baryons and each Higgs doublet, Hip, Qi,2,3, U, D, S, C, B, T, 
with Hi = J0 ft o + fi h - + fi h o + n~ h -), Qi = \{\i UL + fi dL + fiu L + /ijj, U = + fiu R ), 

etc., the interaction terms are 

(r„ + T yA )(-H 2 + Q 3 - T), T w {Hi - Q 3 + T), F hf (Hi + H 2 ), 
T SS (2Q 3 + 2Q 2 + 2Qi - U- D- C- S- B-T). (12) 

Further conservation laws for various unsourced linear combinations give us that Q 2 = Qi 
and U = D = S = C = B. Now putting all the Yukawa interactions in equilibrium gives 
the constraints — H 2 + Q 3 — T = and Hi — Q 3 + T = 0, and using the conservation of total 
baryon number (the weak sphalerons are slow and enter only at the end of the calculation) 
allows the full network to be reduced to just two coupled equations for Q 3 and H = Hi + H 2 : 

-3V h H-3V q Q 3 + T hf H + T ss (28Q 3 -AH) = S H = ~(S Hl + S H2 ) 

+^V q H-9V q Q 3 + T ss (28Q 3 -AH), = 0. (13) 

where T>i = Di-^2+v w -^. We must keep the strong sphaleron and helicity-flip rates explicitly 
because in the limit that either becomes infinitely fast, they drive the left-handed quark 
number, essential for baryogenesis, to zero. The Yukawa interactions, by contrast, simply 
redistribute the asymmetries between the species, transferring the chemical potential from 
Higgsinos, where the source is, into the left-handed quark sector. It is well known that the 
strong sphalerons tend to drive the asymmetry to zero in the massless approximation we are 
working in |TIJ and, since B L = Qi + Q 2 + Q 3 = 14Q 3 — 2H, this is manifest in eq. flT3]). For 
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the helicity flips, which are due to the Dirac mass term \ih\hi that connects the two species 
of Higgsinos, the potential suppression arises from the fact that Sh± = Sh 2 , since hi and 
h\ are the right- and left-handed components of a Dirac fermion, and thus come with equal 
and opposite CP phases, as can be seen in eq. (^). If helicity-flipping interactions were in 
equilibrium this would force Hi = —H 2 and the source for Hi — H 2 would be — Sh 2 = 0, 
leading to the trivial solution for all the particle densities.^ 

The equations fllTf ) can be solved by finding the appropriate Green's function, of the form 
G(z -z') = J G(p)e ip < z - z '^dp/27r, which is a matrix such that (H(z), Q 3 (z)) T = G(z - 
z')(Sh{z'), 0) T dz'. However, we do not need the full matrix G, because we are only interested 
in the linear combination of Q3 and H that gives the chemical potential for left-handed 
baryon number, Bl = lAQs — 2H. The latter can be determined from the single-component 
equation 

/oo 
G B (z-z')S H (z')dz', (14) 
-00 

where G B = — 2Gn + I4G21. Carrying out the thermal averages in the relativistic limit 
(for simplicity of presentation-we do not make this simplification to obtain the numerical 
results), and using the Maxwell-Boltzmann approximation for the distribution functions, the 
source can be written as 

Sh = v w {^j^ ((mV)'(l/£>)' - ^<l/£ 2 )mV) , (15) 

where to a good approximation (pl/E) = T, (1/E) = l/(2T + m), and (1/E 2 ) = (T + 
m) _1 (2T + m) _1 . The appropriate mass m to use here is |-M X ±| (see above eq. @), with 
± equal to sign(/i — m 2 ): this is the appropriate sign for the mass eigenvalue that goes 
continuously into that of the chargino (m = fi) in the symmetric phase outside the bubble. 

The Green's function is determined by the poles of Gb{p), which are the roots of a quartic 
polynomial. For v w < VTD we can solve perturbatively in v w , and the result to 0(v w ) is 

G B (z) =CY,±k ± e- k ^ {l+v w (sign(z)6g ± - 5k ± z)) + 0(v 2 J (16) 
± 

lr This statement is true for the force-term contribution to Sh, but the unequal interactions of Hi and H2 
prevent the spontaneous baryogenesis part of — Su 2 from vanishing. In practice however it is small. 
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where 



kl 



28T ss D h + 3(2T SS + Y hf )D q ± ^j(2%T ss D h + 3(2r ss - T hf )D q f + 16T hf T ss D 2 ) 

lSD q (D h + D q ) 



Sk 3(3D h + AD q )k 2 ± -(3T hf + 3AT ss ) 

± 6D q (6D h + D q )k 2 ± -56T ss D h -Q(T hf + 2T ss )D q { ' 

and C- 1 = 3{6D h + D q ){k 2 + - k 2 _), 5g± = {25k± - l/D q )/k± - 2k±(5k + - 8k-) /{k\ - k 2 _). 
It is important to go to first order in v w here because of a cancellation that takes place at 
0(v ° ): the Green's function at this order is symmetric about the center of the wall, but the 
source is approximately antisymmetric. 

The baryon asymmetry. Having the solution for the density of left-handed baryon 
number, it is simple to find the rate of baryon violation due to weak sphalerons in front of 
the wall: % = — 9T ws n qL = —9T ws T 2 Bl. The time integral of this rate from t = — oo until 
the time the wall passes a given position, where the sphaleron interactions are assumed to 
switch off, can be converted to an integral of Bl(z) in front of the wall, yielding the baryon 
density tib = (9T WS T 2 /v w ) / °° Bi{z)dz. The baryon-to-photon ratio is therefore 

7n B 28350;^/ 



V 



/•OO 

/ B L {z)dz 
Jo 



s 2-K 2 v w g* 

where g* is the number of degrees of freedom at the phase transition temperature. To 



compute / Bl(z)cIz, we use eq. fll4|) and (0), where the z integral can be done analytically, 
leaving a single numerical integral over z'. 

To obtain the final result we must still specify the shape and speed of the bubble wall 
and the values of diffusion constants and interaction rates. For the wall profile, we took 
the kink solution v(z) = (v c /2) (1 — tanh(2:/2u>)) with width w = l/(V\v c ), corresponding 
to the Higgs potential at the critical temperature V(v) = Xv 2 {y — v c ) 2 . We explored a 
range of values of A and v c consistent with the analysis of the phase transition described 
in refs. [O, giving wall widths between w = 0.07 and w = 0.08 GeV -1 , which should be 
compared to the typical critical temperature of T c = 90 GeV. The other parameters varied 
in the ranges 0.003 < A < 0.014 and 120 GeV < v c < 210 GeV. For the diffusion constants 
and decay rates, we take T ss = 6000r„, s |T2|, T ws = a^T (it has recently been argued that 



parametrically r^ s = Coc^T [ 131 , but lattice measurements of the rate are consistent with a 
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Figure 1: Contours of constant 77x0 = 10 10 ?7 in the plane of m2 (the Wino mass) and /1, in 
GeV, u waU = 0.05, 0.1 and 0.2, tan/? = 2, sin 5 = 1 and wall width 2w = U/T. 



value of C ~ l/a w [|rjj so this does not effect the numerical value of our estimate) and we 
have made the rough estimates Thf = 3i/ 2 T/16tc, Dh = 20/T and D q = 3/T. Our results are 
not extremely sensitive to these values. For the wall velocity we took v w = 0.01 — 0.2. 

In general, we find that the baryon-to-photon ratio 77 is sufficiently large only if the CP- 
violating phase 5 of the /1 parameter is nearly maximal, 5 = it/2. In figure 1 we show the 
dependence of 7710 = i] x 10 10 on the chargino mass parameters m 2 and fi for tan (3 = 2 
(favored by studies of the phase transition) A = 0.0065, v c = 160 GeV, and w = 7/T c , and 
several values of the wall velocity. One finds certain parameter values, like m 2 ~ /i when 
v w £ 0.1, where a smaller angle of 5 ~ 0.1 could be tolerated. It is usually assumed that such 
a large phase is excluded by the experimental limits on the neutron electric dipole moment, 
but this can be evaded if the top and bottom squarks are much lighter than their first 
generation counterparts, since these determine the size of the loop diagrams that generate 
the up and down quark EDM's that feed into that of the neutron. This is in some sense 
natural because a light top squark is also needed for a strong enough phase transition to 
satisfy the sphaleron bound. 



Discussion. We conclude by pointing out some of the important differences between 
previous related work and ours, and why we consider the present analysis to be more satisfac- 
tory. The comparison can be divided into two parts: the derivation of the sourced diffusion 
equations, and their subsequent solution and estimation of the baryon asymmetry. With 
respect to the first point we have noted that in the present treatment the source terms and 
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diffusion equations are consistently derived together. In several recent treatments [II], |3|, |J. 
quantum mechanical contributions to various CP-odd currents induced by the wall are de- 
rived, and then converted into sources for diffusion equations using an ad hoc prescription 
introduced in M. If such a prescription could be rigorously derived, the derivation should 
start with a set of transport equations more general than those used here, incorporating 
quantum mechanical corrections. However any such consistent treatment should reduce to 
our equations in the classical limit, and this is not the case, as can be seen from the fact that 
the source terms of |], 0, U] are parametrically different from ours. Furthermore, insofar as 
the thick wall and on-shell quasi-particle approximations we have used are valid, such quan- 
tum mechanically derived source terms should be higher order corrections to those derived 
here. 

One example of the different parametric dependence of our source term concerns the 
variation of tan f3 in the wall. The source terms discussed in refs. Jl|, [3], |^, |6| are proportional 
to v if 2 — v[v2, which is zero if t>i/t>2 is constant in the bubble wall. As pointed out in ref. || 
and confirmed by the methods of ref. | IT| , V1/V2 varies at most by only a few percent over 
the wall in the MSSM. In ref. H, on the other hand, it is simply assumed that t>i/t>2 changes 
by 100% in the wall, which gives an unjustifiably large estimate of the baryon asymmetry. 
We have concentrated on the chargino source because, in our derivation, it is apparently the 
sole exception to this rule in the MSSM. It escapes the suppression because the chargino is a 
Dirac particle, whose mass eigenvalues have a spatially varying phase when Vi/v 2 is constant, 
even though the eigenvalues of M^M have no such phases.0 

There are several other differences between our treatment and previous ones. First, our 

computation has assumed that all the squarks are light and hence the strong sphaleron 

suppression must be taken into account, whereas || |J took all squarks except for t i B and 

b L to be decoupled. In this case the strong sphaleron suppression is evaded 0. Secondly we 

find that for v w > 0.01 (so that electroweak sphalerons are still out of equilibrium on the 

wall passage time scale), the baryon asymmetry is 0(1) + O(10v w ) rather than 0(1 /v w ) as 

2 Actually the suppression would still occur if Winos and Higgsinos had identical interactions, since they 
come with opposite phases. It is the different interaction rates of the two species that undoes the suppression 
coming from constant V\ jvi . 
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found by refs. [|I|, f|]|] This can be traced to the fact that the Higgsino decay rate is taken 
to be zero in the symmetric phase, but non-zero inside the bubbles. In fact we believe it 
should be the other way around since Higgsinos have much more phase space to decay into 
massless quarks and light stops in the symmetric phase than in the broken phase where the 
top quark is heavy. But regardless of phase space, helicity-flipping scattering processes are 
still fast enough in the symmetric phase to give a significant rate of Higgsino damping. 
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